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Abstract 


The D0 collaboration has performed a measurement of the top quark mass 
mt based on six candidate events for the process tt —> bW~^bW~, where the 
W bosons decay to en or /au. This sample was collected during an exposure of 
the D0 detector to an integrated luminosity of 125 pb“^ of ^/s = 1.8 TeV pp 
collisions. We obtain mt = 168.4 ± 12.3 (stat) ± 3.7 (sys) GeV/c^, consistent 
with the measurement obtained using single-lepton events. Combination of 
the single-lepton and dilepton results yields mt = 172.0 ± 7.5 GeV/cP. 
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The pair production of top quarks has been observed in pp collisions at -y/s = 1.8 TeV 
by the CDF and D0 collaborations |]l|. Since the time of observation, the integrated lu¬ 
minosity has more than doubled (to 125 pb“^) and the D0 experiment has substantially 
improved its techniques for measurement of the top quark mass mt- We previously re¬ 
ported a measurement of using events in which one top quark decayed semileptoni- 
cally and the other decayed hadronically (the ‘T -f jets” mode, where £ = e or p), giving 
rrit = 173.3 ±5.6 (stat) ±6.2 (sys) GeV/c^ [^. This letter reports a hrst measurement of m* 
using events consistent with the ti —> bW^bW~ —> b£'^iA)£~V (“dilepton”) hypothesis. This 
independent measurement is important as a direct test of the hypothesis that the excess of 
events over background in both the £ ± jets and dilepton channels is due to tt production. 

The events used in this analysis were recorded by the D0 detector P], which consists of a 
nonmagnetic tracking system, including a transition radiation detector (TRD), surrounded 
by a hermetic liquid argon/uranium calorimeter, segmented in depth into several electro¬ 
magnetic (EM) and hadronic layers, and an outer toroidal muon spectrometer. Electrons are 
identihed using a likelihood method based on the EM shower shape, track ionization, spatial 
match of the track with the EM shower, and TRD response. Muons are required to have 
reconstructed tracks in the central tracking chamber and in at least one of the spectrometer 
layers outside of the toroid, and to have energy deposition in the calorimeter consistent with 
the passage of a minimum ionizing particle. Jets are reconstructed from calorimeter energy 
clustered within a cone of radius Y^^(Ar7)M-~(A0)^ = 0.5 @]. 

We select dilepton top candidate events according to criteria similar to those used in 
our cross section measurement 0. This selection requires two leptons, with the transverse 
energy Et of each lepton > 15(20) GeV for the ep and pp (ee) channels, with |?7e| < 2.5 
and \rj^\ < 1.7, and two or more jets with Et > 20 GeV and |? 7 | < 2.5. In addition, for 
the ee and ep channels we require signihcant missing transverse energy Iprp to discriminate 
against background sources that have no hnal-state neutrinos, while for the pp analysis we 
reduce the Z boson background by rejecting events for which the probability of a £t 
to the Z ^ pp hypothesis is > 1%. We reject much of the remaining background using 
the quantities Ht = Sjets-^T and iJf. = EIt ± Et^Ci), where all jets with Et > 15 GeV 
and \t]\ < 2.5 enter the sum, and Ci is the leading Et electron. The selection requires 
H^{Ht) > 120(100) GeV for the ee and ep {pp) channels. The ee selection reported in 
Ref. 1^ is extended to include an event that contains, in addition to an electron, one EM 
energy cluster without an associated reconstructed track, but with hits in the layers of the 
central tracking system between the interaction vertex and the cluster. This event also has 
a muon near one of its jets, which is evidence of b quark decay, and further enhances the 
probability that the event is an example of tt production. The signal to background ratio 
for such 6-tagged events in which only one of the EM clusters has an associated track is ~ 
15/1. The final sample therefore consists of three ep events, two ee events, and one pp event, 
with expected backgrounds of 0.21 ± 0.16, 0.47 ± 0.09, and 0.73 ± 0.25 events, respectively 
(background sources are detailed in Ref P). Kinematic details for the observed events can 
be found in Ref. 0. 

We reconstruct the events according to the tt dilepton decay hypothesis. After applying 
the invariant mass constraints m{£iVi) = m{£ 2 i' 2 ) = '^w and m{£iVibi) = m(£ 2 i^ 2 ^ 2 ); the 
system remains underconstrained due to the two undetected neutrinos. We supply the 
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needed additional constraints by assuming values for and for two quantities associated 
with the neutrinos (as discussed below). We then solve for the neutrino momenta up to 
a fourfold ambiguity and assign a weight to each solution to characterize how likely it is 
to occur in tt production for the assumed m* [^. We compute the relative weight as a 
function of assumed mt for 80 < < 280 GeV/ c^, employing two weighting schemes with 

differing sensitivities to top production kinematics, decay distributions, and phase space 
volume consistent with the event topology. 

The hrst weighting scheme is called matrix element weighting (WIWT), and is an ex¬ 
tension of the procedure given in Ref. [Q. Here, we require the sum of the neutrino to 
equal the measured We assign a weight 

Woimt) = A{mt)f{x)f(x)P{Eg^\mt)P{Eg^\mt) 

where f{x) is the CTEQ3M p parton distribution function evaluated at the proton (antipro¬ 
ton) momentum fraction x (x) required by the solution, and P{Eg^\mt) is the probability 
density for the lepton energy in the top quark rest frame. The factor A{mt) normalizes the 
average weight of accepted events to unity, independent of the top quark mass. For each 
mass, we add the weights of all solutions. 

The other weighting scheme, called neutrino weighting (z/WT) |^, steps the assumed rj 
for each neutrino through a range of values at each m^. Each step spans an equal fraction 
of the neutrino rj distribution expected in tt production. At each step a weight is assigned 
based on the extent to which the $rp measured in the event agrees with the sum of the 
neutrino pV’s in the solution. The Gaussian resolution of each component of the Ifij, is 4 
GeV. The weights at all rj values are summed to give Wo{mt) in this method. 

In calculating the consistency of the event kinematics with any given mt, we also account 
for detector resolution for jets, leptons, and Ifirp. This is done by fluctuating all the measured 
energies by their resolutions 100 (5000) times for MG (collider data) events and summing 
the weights obtained for each fluctuation. The Gaussian resolutions for electrons and muons 
are ge/E = 15%/a/E ©3% and ct(1/p) = 0.18(p — 2)/p^ © 0.003, respectively, with E (p) in 
GeV (GeV/c). Jets are smeared by double Gaussians designed to model both the inherent 
energy resolution of the hadronic calorimeter (narrow Gaussian) and the contribution of 
large angle gluon radiation to the resolution (wide Gaussian). The Ifirp is then recomputed 
to reflect the changes in jet and lepton energies, and each component is fluctuated with a 4 
GeV Gaussian. 

For both the AJWT and z/WT methods, the up to fourfold solution ambiguity is handled 
by summing the weights for all solutions. There remains a twofold ambiguity in the pairing 
of jets to leptons in reconstructing the event. Both analyses sum over both pairings. In 
addition, initial- and hnal-state gluon radiation (ISR and FSR) can create extra jets that 
further complicate the hnal state. For events with more than two jets, a weighted sum is 
taken over all possible combinations of the three leading Et jets. If jet i is assumed to 
arise from ISR, then the weight is given by exp[—F7’sin6'j/(25 GeV)] |Q. If jets i and j are 
assumed to arise from the same b quark by FSR and have invariant mass mtj, the weight 
assigned is exp[—mij/(20 GeV/c^)]. In each case, the form of the weight is based on the 
characteristics of gluon radiation and the coefficient is determined empirically to maximize 
the sensitivity. The distributions of W{mt) (which corresponds to Wo{mt) after accounting 
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FIG. 1. W{mt) distributions for the six dilepton candidates, using the AfWT (dashed) and 
i/WT (solid) methods. 


for resolution and jet combinations) for the six candidate events, using the two weighting 
methods, are displayed in Fig. |^. 

Due to combinatorics, and the fact that no account has been taken of background, the 
W{rrit) distribution derived from either the A^WT or i/WT method cannot be considered to 
be a probability density. Consequently, a maximum likelihood £t is performed in which the 
W{rrit) distributions for data are compared to the expectations from signal and background. 


The signal is modeled using herwig while the various background sources are modeled 
by ISAJET B , PYTHIA 


and D0 data (for the instrumental backgrounds) 


For both analyses, the maximum likelihood fit proceeds by normalizing the W(rrit) dis¬ 
tribution for each event to unity, and integrating the fractional weights in five 40 GeV/c^ 
bins in rrit. The hrst four bins form the components of a four-dimensional vector Wt for each 
event i. Using the shape of the W{mt) distributions increases the statistical sensitivity of 
the measurement by ~ 25% over a fit to a single-valued mass estimator for each event. The 
likelihood L{mt,ns,nb) to be maximized is 

L = g(nMn. + %) x R , 

Us + Ub 


where and rib are the fitted signal and background levels, gijib) is a Gaussian constraint 
that rib be consistent with expectations, p^rig + rib) is a Poisson constraint that rig + rib be 


7 
























TABLE I. Summary of mt measurements for full and partial data sets. Uncertainties are 
statistical only. 


Channels Eit 

MWT 

nWT 

eiJ, + ee + pp 

168.1 ± 12.4 GeV/c^ 

169.9 ± 14.8 GeV/c^ 

ep -|- ee 

167.9 ± 12.6 GeV/c^ 

173.2 ± 14.0 GeV/c^ 

ep 

173.1 ± 13.3 GeV/c^ 

170.1 ± 14.5 GeV/c^ 


consistent with the sample size iV, and fs and fb are the fonr-dimensional probability den¬ 
sities for signal and background. The probability densities f{w) are estimated by summing 
the contributions of four-dimensional Gaussian kernels placed at the location of w for each 

iH 


event in the signal MC or background samples 


Using the estimated fg and fb, L is 


maximized with respect to Ug and Ub for each value of rrit at which we have generated MC 
events. The maximum likelihood estimate of mt {mt) and its error (d) are determined by a 
quadratic fit to — In L for the nine points about the minimum. 

Applying the maximum likelihood fit to the data, we determine the top quark mass to 
be mt = 168.1 ± 12.4 GeV/c^ (AfWT), and mt = 169.9 ± 14.8 GeV/c^ (z/WT), where the 
uncertainties are statistical only (see Fig. The results of hts to subsamples of the data 
are listed in Table |. 

To study the properties of our mass estimator, 1000 MC “experiments” are generated 
by randomly selecting three e/x, two ee, and one fifi events from simulated top (with input 
mass mf^^) and background samples, according to the estimated background contamination. 
Due to the small sample size the estimator is not normally distributed. Figure ^ shows a 
typical distribution for m*. The distribution is characterized by its median and width (half 
of the shortest interval that contains 68.3% of the experiments). The width of a Gaussian 
fit to the pull distribution {nit — mf^^)/a is consistent with unity. This verifies that a is 
an unbiased estimate of the statistical uncertainty. The pulls means can differ from zero 
because of the non-Gaussian tails of the frit distribution. Table II lists medians, widths, and 
means and widths of the pulls for different m^^. The properties of frit are very similar for 
both methods. 

The sources of systematic uncertainty are summarized in Table B A major uncertainty 
arises from the jet energy calibration, which proceeds in two steps. The hrst step uses y-l-jets 
events to relate the well-calibrated electromagnetic scale to the hadronic scale. Effects con- 


TABLE II. Properties of the maximum likelihood estimate rhf. 



^MC 

GeV/c^ 

Median 

GeY/cf 

Width 

GeV/c2 

Pull 

Mean 

Pull 

Width 


160 

161.6 

15.8 

0.12 ± 0.03 

1.03 

MWT 

170 

172.2 

16.7 

0.11 ± 0.03 

0.99 


180 

180.5 

17.3 

0.00 ± 0.03 

0.98 


160 

161.5 

14.4 

0.17 ± 0.03 

0.96 

nWT 

170 

172.2 

16.2 

0.08 ± 0.03 

0.98 


180 

180.5 

18.1 

0.03 ± 0.03 

1.03 
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FIG. 2. Sum of the normalized candidate weights grouped into the five bins considered in the 
maximum likelihood fit (circles) for the (a) A^WT and (b) i^WT analyses. The uncertainty on 
these points is taken from the RMS spread of the weights in MC studies. Also shown are the 
average weights from the best-fit background (dashed) and signal plus background (solid). The 
— InL distributions and quadratic fits are inset. 


sidered include the overall hadronic response, energy added to jets by multiple interactions 
and uranium noise, and the spread of the hadronic shower outside of the jet cone H]. The 


second step follows from a detailed comparison of 7 + jets events in MC and collider data, 
providing a correction that depends on jet 7 and ensures that the energy scale in data 
matches that in MC. An additional correction is applied to jets having a muon within the 
jet cone. Under the assumption that the muon was produced in the semileptonic decay of a 
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FIG. 3. Distribution of fht for “experiments” with = 170 GeV/c^ (marked by vertical 

line) using the A^WT (dashed) and i^WT (solid) analysis methods. 


TABLE III. Systematic errors in the measurement of mt. 

Source 

Jet energy scale 
Signal model 
Multiple interactions 
Background model 
Likelihood fit 
Total 


Error (GeV/c^) 
2A 
1.8 
1.3 
1.2 
1.3 
3T 


heavy quark, the jet energy is increased take into account the muon and neutrino energies. 
After these corrections, we hnd (using the y+jets and smaller 2’(—>• ee) + jets samples) that 
the data and MC scales agree, with an uncertainty of 5{Et) = 0.025i?T + 0.5 GeV. MC tests 
are run on samples with the jet energies rescaled by ±5 to yield an energy scale uncertainty 
of 2.4 GeV/c^ in m*. 

Other sources of uncertainty arise from differences among models of tt and background 
production, multiple interactions, and the likelihood £t procedure, isajet is used as a cross¬ 
check of the HERWIG tt production model. The effect of multiple interactions is estimated 
using MG samples with additional random pp interactions overlaid on tt production. The 
contributions from all sources are summed in quadrature to give the systematic uncertainty 
on the measurement (see Table B 

Taking account of the 77% correlation between the AJWT and i^WT analyses, we mea¬ 
sure the mass of the top quark in the dilepton channel by combining the two results: 

rrit = 168.4 ± 12.3 (stat) ± 3.7 (sys) GeV/c^. 

This value is in good agreement with the measurement from the i + jets channel p|, con¬ 
sistent with the tt hypothesis for both channels. We combine the results of the single- and 
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dilepton analyses by propagating the systematic uncertainties in each channel with corre¬ 
lation coefficients of either zero (for MC statistics and background model) or unity (for jet 
energy scale, multiple interactions, and tt production models). Using all dilepton and i + 
jets tt candidates, we measure: 

rrit = 172.0 ±5.1 (stat) ± 5.5 (sys) GeV/c^, 

or, combining statistical and systematic errors in quadrature, = 172.0 ± 7.5 GeV/c^. 
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